Analysis of mechanical performances of cylindrical wall in 3D printing processes 
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In order to enhance the overall structural stability of 3D printed concrete wall, we propose a novel scheme to 
produce a wavy wall with curvature along its contour direction. To validate the idea of the scheme, a single wave 
wall is set as a cylindrical wall and its mechanical performances are analyzed. The mathematical 3D printing 
model of the wavy wall is formulated by the shell theory while taking into account of additional parameters of 
the printing processes, the model will be used to analyze the two failure mechanisms of the cylindrical wall: 
elastic buckling and plastic collapse. Compared with the results of Suiker’s straight wall, it is found that when 
the parameters were same, the stability of the cylindrical wall is more than twice of the rectangular wall. Our 
studies indicate that it’s a feasible scheme to improve the printed structural stability by increasing curvature. 
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I. INTRODUCTION 


The conception of 3D printing, also referred to as addi- 
tive manufacturing (AM), rapid prototyping (RP), or solid- 
freeform technology (SFF), was developed by Charles Hull. 
Modern additive manufacturing technology was introduced 
approximately three decades ago, this technology has devel- 
oped particularly rapidly and achieved great success. The ap- 
plications include biomedical engineering, automotive engi- 
neering, mechanical engineering, aerospace engineering, food 
engineering, building civil engineering, etc. [1-7]. The term 
” Additive Manufacturing” was finally chosen by ASTM to 
distinguish it from traditional ’Subtractive Manufacturing”. 
The main form of 3D printing is to establish a 3D model digi- 
tal file, use software to slice the model (the accuracy depends 
on the printing resolution), and finally use the printing equip- 
ment to print layer by layer [8]. The advantages of AM tech- 
nology over traditional manufacturing are: it can manufacture 
parts with complex geometric shapes, reduce material waste, 
reduce energy consumption, increase customer flexibility, etc. 
[9, 10]. 


In 2009, Bourell et al. [11] published a roadmap for AM 
based on a workshop of 65 key people in AM. Their report 
explored important facets of the AM including: (1)design, 
(2)process modeling and control, (3)material, processes and 
machines, (4)biomedical applications, (5)energy and stain- 
ability application. As for the diversified development of addi- 
tive manufacturing (AM), Behrokh Khoshnevis proposed the 
*contour crafting”, also known as 3D concrete printing tech- 
nology. Due to the limitations of materials and economic fac- 
tors, the real potential of this technology has not been ex- 
plored. However, in some extreme conditions, it effectively 
compensates for the weaknesses of traditional buildings, such 
as the application of mobile cabin hospitals and the construc- 
tion of habitats on other planets [12-14]. At present, the anal- 
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ysis of mechanical performances in the 3D printing processes 
is relatively rare, and the main research focuses on the mate- 
rial science [15-18]. Due to the complexity and diversity of 
parameters involved in the 3D printing processes, including 
the geometric features of printed objects, the impact of time 
effect on material strength, the inhomogeneity of self-weight 
load, the moving speed of printing nozzle, and the inevitable 
geometric defects, the printing parameters are usually set in a 
trial-and-error manner. There are few researches on the me- 
chanical performances of concrete in the printing process. 


According to the research of Suiker et al.[19—23], elas- 
tic buckling and plastic collapse are the main mechanism of 
structural failure in 3D concrete printing processes. The con- 
tribution of Suiker is to analyze the mechanical performances 
in the 3D printing processes of straight wall structures, and 
give the corresponding mechanical model for the first time. 
It was found that the failure behavior of the structure in the 
printing process could be analyzed well, the number of print- 
ing layers of the straight wall structure could be predicted, 
and the theoretical guidance was given to find the best set of 
printing parameters. The results show that the stability of the 
rectangular wall structure is still insufficient, that is, the struc- 
ture buckling rather than collapse when the plastic limit of 
the material is reached. Stability studies are of great impor- 
tance in improving material utilization, rapid prototyping and 
strengthening structural bearing capacity. In particular, how 
to improve the overall stability of the printed wall structure is 
a key scientific problem when the printing speed is constant 
and 3D printing materials such as concrete, which have strong 
fluidity and slow curing, are used. 


In order to solve this problem, a simple and feasible uni- 
versal scheme is proposed in this paper, that is, the curvature 
of the straight wall along its contour direction is appropriately 
increased to print a wavy shape, the scheme is illustrated in 
Fig.l. Each ripple of wavy wall, its can be cylindrical or not, 
its curvature also is variable. In order to verify our scheme, 
this paper studies the mechanical performances of cylindrical 
wall during 3D printing processes [24-26]. 
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Figure 1: Wavy wall schematic diagram. 


The manuscript is organised as follows. In section 2, the 
governing equations of buckling of the cylindrical wall are 
derived. In section 3, two possible failure mechanisms of 
cylindrical wall in the 3D printing processes were described. 
In section 4, the numerical calculation method was used to 
obtain the buckling results under the influence of two curing 
parameters, and the Suiker number was used to evaluate the 
competitive relationship between the two failure mechanisms 
of the cylindrical wall. In section 5, the same printing pa- 
rameters as Suiker [19] were used to compare and analyze the 
structural stability of cylindrical wall controlled by curvature. 


I. HETEROGENEOUS CYLINDRICAL WALL 
SUBJECTED TO SELF-WEIGHT 


A cell of wavy wall is taken for analysis, whose curvature 
can be changed, that is, k = k(r), r = variable. However, 
when the curvature is constant, i.e. r = constant, and the 
structure is closed, an cylindrical wall with maximum struc- 
tural stability is obviously obtained, see Fig. 2. 
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Figure 2: Cylindrical wall evolution: a straight wall(left), a wavy 
wall(middle), a cylindrical wall(right). 


The mechanical performances of cylindrical wall will be 
heterogeneous in the vertical direction due to the influence 
of material curing behavior in the 3D printing processes. In 
addition, the cylindrical wall subjected to self-weight, and the 
stress in the vertical direction (self-weight direction) is also 
non-uniform [27-29]. Therefore, based on the equilibrium 
equation and boundary conditions, the mechanical model of 
the heterogeneous cylindrical wall under the action of non- 
uniform in-plane forces should be described. In section 2, the 
differential equation of buckling is derived by the principle of 
potential energy minimization. 


A. Energy formulation 


Figure 3: Cylindrical wall under self-weight. 


Cylindrical wall radius is R, height is L, uniform thick- 
ness is h, see Fig.3. In axisymmetric condition, u(x, 0) = 
u(x), v = 0, w(x,0) = w(x). According to the shell theory 
[? ? ], the strain-displacement relationship can be expressed 


as: 
w 
Ex = Ug, Ey = R Eay = 0, 
Xa = W rz; Xy = 0, Xaey = 0. (1) 
in which us = 2 = ËW, yis the longitudinal 
in which us = 5%, Wrs = Gæ; «x is the longitudina 


coordinate; 0 is the circumferential coordinate; Ez, Ey, Exy 
are the linear strain components of the shell middle surface; 
Xe» Xy, Xey are the curvature changes; and u, v, w are the 
displacements in the longitudinal, tangential and radial direc- 
tions, respectively. 

The total strain energy of shell U is composed of the mem- 
brane strain energy U, and the bending strain energy Up. 


U =U; +U, (2) 


where 


Eh 2, 2 1 
s= f + QE x 
U. || z E Ey + 2VEzEy + 


—v 
=y Ewy dydz, 


(3) 


D, 
m= f | FE DE + + wxx +20 - x2] dye, 
(4) 


Here, the asterisk subindex , indicates that the material pa- 
rameters are heterogeneous in space. F, is the non-uniform 
stiffness modulus, the heterogeneous bending stiffness is 
given by 


Eh’ z 


D= D, (x,y) = ia — v2) E, = E, (x,y). (5) 


The work done by the self-weight during buckling of wall 
can be expressed as 


w=3 f fone- (E + Ey] uae, ©) 
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Combining Eqs.(2) and (6), the total potential energy II of 
the shell follows as 


[e? +e? 4 
ES ae a [ez + Ey T 


a —— e2 ]dydz 


sfe [xa + Xo + WXeXy 


+2(1- v)Xayldydz 


/ J: 5 pgh(x 


QExEy 


l) [v?, +u? p] dydz. (7) 


B. Equilibrium equation and boundary conditions 


The equilibrium equation of an inhomogeneous cylindrical 
wall subjected to self-weight can be obtained by minimizing 
potential energy by Eq.(7) 


1 
= 3 | [> [Wr OW rx + W yyOW yy + VW x2 OW yy 


2(1 — v) 
R2 


E, ia 1 v 
+3] e xÔU s + p wow + pe’ 


7 —v) 
R2 


1 
=5 I fia — 1)v,260,2 + W,2dW,2|dydz, 


+ VW yyOW re + W xyðW xy |dydx 


+ pcdw + UyOuy|dydx 


(8) 


Inserting Eq.(1) into Eq.(8) and the displacement variational 
principle is adopted when dealing with axisymmetric buckling 
of cylindrical wall. The differential equation follow as 


1 


E (pgh(a a I)Wwa) o = 


(9) 


In addition, the boundary conditions in the edges x = 0 and 
x = lare 


wļ|s=0 = 0, W,2|2=0 = 0, 


ee ee = 0, (DiW tu) alei = 0. (10) 


Til. FAILURE MECHANISM OF CYLINDRICAL WALL IN 
3D PRINTING PROCESSES 


plastic collapse 


J 


elastic buckling 


Figure 4: Failure mechanism of cylindrical wall in 3D printing pro- 
cesses. 


Two failure mechanisms may occur during 3D printing: 
elastic buckling and plastic failure, see Fig.4. When the ma- 
terial hardening rate is slow, large deformation will occur at 
the bottom of the structure. According to the experimental 
[? ], the first layer will be less deformed under constraints, 
while the final thickness of the second layer will reach more 
than twice that of the first layer when plastic collapse, which 
can be called “Elephant foot deformation”. The reasonable 
utilization of this phenomenon can enhance the structural sta- 
bility to some extent. 

In section 3, the governing equations of these two fail- 
ure mechanisms are established, and the competitive criteria 
for elastic buckling and plastic failure are given by introduc- 
ing the time effect factors related to the curing processes and 
printing. 


A. Time effects on buckling of cylindrical wall 


Figure 5: 3D printing processes of cylindrical wall. The bottom is 
the Lagrangian coordinate system x, the print head is the Eulerian 
coordinate system X. 


As the cylindrical wall length L grows with the printing 
process, this processes is simplified as the growth of whole 
cylindrical wall structure along x—direction at a constant 


speed 4 dt — |, Where the speed l can be obtained from the 
flow donb Q= lu, hTh. 
— (11) 


Un hT, 


In which Q is the volume of printing material discharged from 
the nozzle per unit of time, vn is the horizontal moving speed 
of printing nozzle, T; is the time required to print a layer, h 
is the thickness of the shell. The printing processes can be 
determined by using the Eulerian coordinate system attached 
to the end of the printing nozzle, see Fig.5. The relationship 
between Eulerian coordinate X and Lagrangian coordinate x 
is l = x + X. The Eulerian coordinate X is as follows 


X = X(a,t)=2-1, (12) 


In general, the printing speed can be considered constant, that 
is, 1 = constant. And set the time as /(0) = 0, and it matters 
l = lt. Eq.12 can be simplified as 


X = X(a,t) =2— it, (13) 


Where t represents time. If the printing speed is variable the 
problem is more complex and will be considered in future 
studies. 

Consider the material at point x = 0, the evolution of elas- 
tic stiffness during curing can be expressed as 


E, (x = 0,t) = g«(t) Eo, (14) 


In which g, = g,(t) is curing characteristic function, Eo 
is the initial stiffness of material. As shown in Suiker [? ], 
two most representative curing functions in Eulerian coordi- 
nate system can be obtained. 

(a)Linear curing function: 


(X) =1-X, (15) 
(b)Exponentially-decaying curing function: 


# (X) <1 + 1 —y)exr(®), ae 


Where 
X 
, ên € {éh Eb} (17) 


The bar on the parameter indicates that it is a dimensionless 
parameter. 


Combineing Eqs.13 and 17, X (x, t) follow as 


X =X (x,t) = Ha- Ept, 


x € {£k ép} (18) 


Substituting Eq.(18) into Eq.(9), the dimensionless buckling 
equation can be obtained. 


— a2. 121 =") ha l__ 
(Jw XX) XX ra X4 -3 (=) g JW 
r Ep h eee 
i 
=æ (XUx) x = 0, 
EE 
(19) 


The dimensionless parameter W, T, h, ler and E pare 


w= w/h, 

Do e r= (afii, i= (afin, 
Tor = Enr = (Se, w= SE 

En = (oye, gne (ehti) 


(20) 


Where a is the dimensionless amplification factor. The pur- 
pose is to change a@ to make h = 1 when the parameters are 
given. Eq.(19) can be transformed into 


o k 
(G0 xx) xx tV X +121- v’) g9 T 
cE Se 
S 
23 (XW) x =0, on 


Where k is dimensionless curvature, the buckling equation 
will be controlled by curvature k and curing rate €. 


= 1 
k=-. (22) 
T 
Inserting Eq.(18) into Eq.(10), the boundary conditions 
along X = —« and X = 0 are as follows: 
rel, oe = 0, Uz|X=-x = 0, 


B. The solution of eigenvalue problem 
Eq.(21) can be solved by Galerkin method. 
0 
ie (Rôwn)dX = 0, (24) 
X=—K 


with the residual R as 


R= R(X) = (9.0 xx) xx tvs X 
fn 
F 1 
+ 12(1 — v?)—--9,w — = (Xx) x =0 (25) 
ÉE EE 


In which ôw, is test function, the solution W, = 
5 Cox and dW, = (0W/OC;,,)dC,. Cn repre- 
senting generalised coordinates. Eqs.23 and 24 were used 
to obtain generalized coordinate algebraic equations, and 
Newton—Raphson iteration method was used to solve the 
equations numerically. In essence, it describes the eigenvalue 
problem of axisymmetric buckling of cylindrical wall. 
Consider the influences of dimensionless parameters in 
Eqs.(20) and (22), i.e, diameter-thick ratio r/h, the thickness 
h, buckling length ler, curing rate, E g, amplification factor 
a, vertical printing speed i. The generalized buckling length 


function can be obtained, ler = ler (Ep, K). 


C. Plastic collapse of cylindrical wall 


The cylindrical wall will be failured when the plastic yield 
strength o, of the material is reached. The subscript p stands 
for “plastic collapse ”. Generally speaking, the point location 
of plastic collapse is at the bottom of the cylindrical wall, that 
is, at x = 0, where the axial stress caused by self-weight 
reaches its maximum value. The yield strength can usually be 
expressed as 


pgl = |p|, (26) 


Where |.| on behalf of the absolute value of the yield strength 
Op, lis the length of the cylindrical wall. It is treated similarly 
to the change of bending stiffness D,. The yield strength at 
the bottom varies with time can be set as 


p(x = 0,t) = hy (top 0. (27) 


Where cp o is the yield strength at the moment when t = 0 
in the printing processes, h, = he (t) is the curing character- 
istic function of yield strength evolution over time. Since the 
plastic yield limit is determined by the material and the curva- 
ture of the geometric parameter has no influence on it, dimen- 
sionless plastic collapse length changes with two type curing 
functions which are exactly the same as those in Suiker [? ]. 
(a) Linear curing yield function: 


= 1 =] 
lp, = = 
I= Eg 


<1, (28) 


(b) Exponentially-decaying curing yield function: 


lp — [Yo + (1 — yo) exp(—E5lp)] = 0, (29) 


Newton-Raphson iterative method was used to solve Eq.(29), 
and a closed numerical solution with accurate approximation 
was obtained. where the R?-value of this approximation lies 
between 0.993 and 0.999 [? ]. 


a 5-1 
lp = lpo | 1+ ae , 
i 
= & 1.181 
rar = Eref Ye) = 140.8444,’ 
P = B(yo) = 1.466(y,)°"””, (30) 


Plastic yield length l; and curing rate €, can be expressed as 


L= pglp 
P |op,ol’ 


Z, _ §0|2p.0l 


Es E {ELES}. (31) 
pgl 
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Figure 6: Curve of plastic yield length 1, with curing rate € = E 
under the influence of linear curing parameters. 
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Figure 7: Curve of plastic yield length lp with curing rate € = E 
under the influence of exp-decaying curing parameters. 


è : : : zl 
As shown in Eq.(31), dimensionless curing rates €, can 
be calculated priori from parameter of material and printing 


process. Such as Eq.(28), when E — 1, yield length l close 
to infinity. Therefore, for linear curing processes, the value of 
curing rate approaching 1 should be considered as the upper 


limit of plastic collapse. When E€ > 1, the yield strength 
increased faster than the stress controlled by the printing speed 
due to the curing rate €4. Therefore, plastic collapse will not 
occur. 


As shown in Eq.(30), when E — œ, the plastic collapese 
length ly converges to Ip,0°7- Therefore, for the exp-decaying 
curing processes, Lo - y should be considered as the upper 
limit of plastic collapse. 


D. Suiker number and Suiker criterion for elastic buckling 
and plastic collapse 


Plastic collapse occurs when the plastic yield length is less 
than the elastic buckling length lp < ler. Conversely, when 
lp > ler, elastic buckling is the failure mechanism. This stan- 
dard for analyzing possible failure mechanisms can be conve- 
niently represented by calling geometry, material, and printing 
parameters, and can be obtained by using Eqs.(20) and (31) 


= <A elastic buckling, 
p 
l er ea . 
A >A plastic collapse, (32) 
p 
Where 
x h_\1 |op,0l 
A=(a—)3—*), 
Do” (pg)s 
l= lela tele): (33) 


The dimensionless parameter A mentioned in the Eq.(33) 
was first proposed by Suiker [? ], because of its importance, 
A was called Suiker number in this paper. A is constant when 
the geometric parameters of the printed object and material 
parameters are given, and does not change with time. When 
evaluating the competitive relationship between elastic buck- 
ling and plastic collapse of straight wall structures, A can be 
calculated in advance so as to predict the failure mechanism. 
In printing process, the method of trial and error is often used 
to analyze the failure mechanism. Even under the same ma- 
terial parameters, the failure of structure will also change due 
to the difference of geometric parameters. The Suiker number 
A provides a clear direction to achieve the maximum mate- 
rial utilization by improving material properties or changing 
geometric parameters, which is of great significance for an- 
alyzing the mechanical performances and failure mechanism 
of straight wall structures in 3D printing processes. 


IV. NUMERICAL SOLUTION OF CYLINDRICAL WALL 
FAILURE PROBLEM 


In the result analysis, the critical buckling length Ter (€p, k) 
can be obtained from Eqs.(19) and (23). Substituting wW into 
Eq.(24) and transform it into an eigenvalue problem. The min- 
imum eigenvalue is obtained to represent the critical buckling 
length Ter and is calculated by assuming that the displacement 
is polynomial. 


w(X) =) OX n=1,2, N. 6A 


Where Cn represents unknown generalized coordinates. By 
using the Suiker number given by Eq.(33), the competitive 


relationship between the critical buckling length ler(Ep, k) 


6 


and the plastic collapse length E) in the linear curing pro- 
cesses and the exp-decaying curing processes can be evalu- 
ated. 


A. Numerical solution analysis under the influence of curing 
processes 


When h = 1, the buckling length depends on two dimen- 
sionless parameters, the curing rate €p and the curvature k, 


eS a tee, tet 
i.e., ler(Ep, k). Poisson’s ratio is v = 0.3, and set the range 
of curvature is J <k< F, 


— k=1/5 
309 : p ~ 1⁄6 
linear curing = 1/7 
— 18 
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——. 1/10 
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(a)Effect of linear curing rate 
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(b)Effect of exp-decaying curing rate 


Figure 8: In the 3D printing process, buckling length Ter changes 


yiti curing rate € = Ep, where the range of curvature is Z <k< 


5° 


For a cylindrical wall, Fig.8 illustrates that the buckling 
length ler increases with the increase of curing rate € when 
the curvature k is given; the buckling length ler with curva- 
ture k increases when the curing rate £ is given; the buck- 
ling length Ter has a larger increase amplitude when the cur- 


vature k is greater; the buckling behavior of homogeneous 
cylindrical shell under self-weight is described when curing 


rate Es — 0. The difference between the two curing meth- 
ods is that the linear curing grows faster and the exp-decaying 
curing converges to lero - y with E > o. 

This describe the variation buckling characteristics of cylin- 
drical wall subjected to self-weight from homogeneous to het- 
erogeneous during the printing processes, as well as the influ- 
ence of curvature on structural stability, and conform to the 
general fact that buckling length should increase with the cur- 
ing processes. 


B. Numerical analysis of elastic buckling and plastic failure 


This section mainly analyzes the competitive relationship 
between elastic buckling and plastic collapse. 

Under the condition of linear curing process, curvature 
k= z is selected, and the numerical solution with accu- 
rate approximation was obtained. Where the R?—value of this 


approximation is 0.999. 


Ten = 8.7192 + 3.0652 exp(18.2639E 9), (35) 


When és — 0, that is, the buckling length is independent of 


curing rate, lero = 11.65, the competitive relationship can be 
described by Eqs.(28) and (35). 


|g" 


e _ (8.7192 + 3.0652 exp(18.2639€75))(1 — BE), (36) 
P 
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Where 8 = E, /Ẹp. 
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Figure 9: The failure mechanism of cylindrical wall during the linear 
curing processes. It is the curve of length scale ratio changing with 


curing rate € = ES under selected stiffness curing ratio. Where the 
range of 6 = E JE from 5 to 10. This curve is plot by Eq.(36), A 
represents the competitive relationship between elastic buckling and 
plastic collapse during linear curing. 


The values of 8 are 5,6,7,8,9,10, see Fig.9. Specifically, 
plastic collapse occurs when the point on the curve in the fail- 
ure mechanism graph is greater than the value of A, and oth- 
erwise elastic buckling occurs. The competitive relationship 
obviously depends on the change of 8. As ( gets bigger, the 
length scale ratio ler / l; drops faster. For Fig.9, when 8 > 6, 
as long as A> ler,0 = 11.65, plastic collapse will not occur 
and elastic buckling is the only failure mechanism. 

When A < Lea = 11.65, plastic collapse will occur when 


the curing rate Ei is small, and elastic buckling will occur 
when the curing rate a is large. As Er increases, ler / la de- 
creases along the negative axis of ler / ls until ler / l; = 0, and 
elastic buckling is the only failure mechanism. Specifically, 
elastic buckling is the only failure mechanism when 8 =5 
and A > (ler/lp)max = 14.58. When 0 < A < (ler /lp)maz, 
the range of 0 < E < a,b < Ep is elastic buckling, the 
range of a < č < bis plastic collapse, a, b can be solved 
abo o SA — 
by Epg(ler/lp), where (ler/lp) = A. 
Under the condition of exp-decaying curing process the nu- 


merical solution with accurate approximation was obtained. 
Where the R?—value of this approximation is 0.999. 


ler = 22.9168 — 11.4916 exp(—8.0559€). G7) 


Combined with Eqs.(30) and (37), the competitive relation- 
ship can be expressed as 


[S 


r (22.9168 — 11.4916 exp(—8.0559E) ) 
p 


~l 


—1 


1+ ———, . (38) 
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Figure 10: The failure mechanism of cylindrical wall during the exp- 
decaying curing processes. It is the curve of length scale ratio chang- 
ing with curing rate € = & A under selected curing ratio. Where the 


range of 8 = E ie from 0.1 to 40. This curve is plot by Eq.(38), A 
represents the competitive relationship between elastic buckling and 
plastic collapse during exp-decaying curing. 


The values of 8 are 0.1,0.2,0.5,1,2,5,10,20,40, see Fig.10. 
The specific analysis is similar to the above. When 8 > 20, 
as long as A > 11.65, plastic failure will not occur, elastic 


buckling is the only failure mechanism. When A < 11.65, 
the range of a < én < b is elastic buckling, the range of 
0< ep <a,b< és is plastic collapse, a, b can be solved by 


al o o o 
€p(ler/lp).- a O 

In particular, when 6 = 0.1 and A > (ler/lp)max = 
22.54, elastic buckling is the only failure mechanism. When 


11.65 < A < (ler/Lp)max, the range of 0 < és <a,b< a 


is elastic buckling, the range of a < En < b is elastic col- 
lapse. 

It can be seen from the analysis in Fig.9 and Fig.10 that 
the stiffness curing ratio ( has a great influence on the com- 
petitive relationship between elastic buckling and plastic col- 
lapse. When £ is smaller, it is more prone to plastic collapse; 
On the contrary, elastic buckling is more likely to occur. The 
intersection of the length scale ratio curve and A represents 
the critical position of the two failure mechanisms. Accord- 
ing to the change trend of the curve, during the linear curing 
process, when 8 = 5, the length scale ratio ler/lp first in- 
creases and then decreases; And when 8 > 5, lox /lp is de- 
creasing monotonically. In the exp-decaying curing process, 
when 3 < 5, ler j l first increases and then decreases; When 
B > 10, lo / b first decreases and then increases. Thus, it can 
be seen that the length scale ratio 6 has a great influence on 
the change trend of the curve, that is, under the given geomet- 
ric parameters, ( determines the dominant position of the two 
failure mechanisms of elastic buckling and plastic collapse. 


parameter value 
Thickness of wall h = 43.5(mm) 
Height of layer tı = 9.2(mm) 


p = 2020(kg/m*) 
Un = 83.3(mm/s) 
v = 0.3(—) 


Density of material 
Printing nozzle speed 
Poisson’s ratio 
Linear curing 
Initial stiffness modulus Æo = 78100(Pa) 

Curing rate a = 0.02, 0.05, 0.1(—) 
Exp-decaying curing 
Initial stiffness modulus Æo = 78100(Pa) 
Eig = 0.02, 0.1, 0.2(—) 


ye = Fe =2(-) 


Curing rate 
Curing stiffness ratio 


Table I: Printing curing parameters of Suiker 


V. INFLUENCE ANALYSIS UNDER GIVEN PRINTING 
PARAMETERS 


In this section, the buckling behavior of cylindrical wall 
has been simulated by commercial finite element software 
ABAQUS during the 3D printing processes, and the above 
model was verified. Then, the same printing parameters as 
Suiker [? ] were used to compare the cylindrical wall with the 
b = 625mmd = 250mm rectangular wall. Finally, the fail- 
ure mechanism, the competitive relationship between elastic 


buckling and plastic collapse, of the cylindrical wall printed 
by this kind of concrete is analyzed. 


A. Elastic buckling 


=> => => 
initial growth buckling failure 


Figure 11: Buckling growth schematic of cylindrical wall. 


Fig.11 depicts a schematic diagram in which the cylindri- 
cal wall grows at a constant rate in the direction of X axis 
and the structure becomes unstable from the initial small de- 
formation to the critical buckling height until the structure 
fails. The parameters of cylindrical wall are listed in Ta- 
ble I, and the mechanical behavior of buckling under self- 
weight is analyzed. Two common curing characteristics, lin- 
ear curing and exponential attenuation curing, are selected 
for simulation. According to the Eqs.(14) and (20), č = 
0.02, 0.05, 0.1 was selected for the linear curing process and 
Ea = 0.02, 0.1, 0.2 for the exp-decaying curing processes, 
so as to verify the parameter model. 


c numerical solution 


307 = — - - FEM solution 
linear curing k=1/6.4 


buckling length L. 


0.0 0.2 0.4 0.6 0.8 1.0 
displacement w 


Figure 12: Results comparison diagram between numerical and FEM 
under the influence of linear curing parameters. 


exp-decaying curing k=1/6. 4 numerical solution 


— - - FEM solution 


buckling length 7,, 


T T 
0.0 0.4 0.8 L2 1.6 2.0 
displacement w 


Figure 13: Results comparison diagram between numerical and FEM 
under the influence of exp-decaying curing parameters. 


Note that, in order to fit the actual printing size and fa- 
cilitate numerical calculation, œ in Eq.(20) is amplification 
factor. When the actual thickness of wall is 43.5mm, a = 
8.297 to ensure the dimensionless thickness of wall h = 1. 
Where, under the given parameters, the dimensionless curing 
rate € and curing rate € are positive linear relationship. 
By Eq.(33)4, when dimensionless thickness of wall h = 1, 
dimensionless buckling length ler is determined by the di- 


mensionless curing rate Ep and dimensionless curvature k, 
namely ler = Le (Ep, k). 

The results of bifurcation buckling obtained by FEM are 
relatively close to parameter model, see Fig.1213. In partic- 
ular, under the influence of linear curing, when E = 0.02, 
the difference is 3.1%; when č = 0.05, the difference is 
8.6%; when E = 0.1, the difference is 25%. Under the 
influence of exp-decaying curing, when Ee = 0.02, the dif- 
ference is 1.9%; when č = 0.1, the difference is 7.2%; 


when Er = 0.2, the difference is 10%. 

When width of rectangular wall is b = 625mm,d = 
250mm, curvature should be k = 2rd) œ 1/6.4, 
where, the thickness is 43.5mm, the layer height is 9.2mm. 
According to Eqs.(11)(14) and (20)4, can get € = 
0.000256(s~),i = 0.000438(m/s) and i. = 0.0254. For 
this concrete, when the range of curvature is ł <k< b 
the range of dimensionless curing rate increased from 0.0199 
to 0.0397. Within this range, the result of bifurcation buck- 
ling simulated by FEM is very close to that of the parametric 
model, which indicates that the parametric model can simu- 
late the bifurcation buckling behavior of 3D printing concrete 
cylindrical wall well. 

In the classical theory, the solution of the governing equa- 
tion for the buckling of cylindrical wall is written in the form 
of series approximation by trigonometric function: w = 
X ` A- sin(ax)-sin(GBy), and the solution is also sinusoidal 
along the circumference, which is more in line with the actual 
situation. In the model derivation in section 1, considering 


that the initial stiffness modulus Ep is low, it is impossible 
for the 3D concrete printed cylindrical wall to remain intact 
under large deflection, the parameter model is constructed un- 
der the axial symmetry condition, and the solution is set as 
w = w(x). When the value of ga is large, the cylindrical 
wall has a higher buckling length, the greater difference in its 
horizontal buckling form, and the significant difference in the 
results. 


B. Comparison and analysis 


According to b = 625mm, D = 250mm rectangular wall 
of Suiker [? ], the upper bound of the failure mechanism 
is determined by plastic collapse, the lower bound is deter- 
mined by the elastic buckling, where the parameters are listed 
in Tablel, the curing characteristics is determined by equa- 
tion (39)(40). When used with identical parameters of Suiker, 
radius is r = 278.4mm , thickness is h = 43.5mm and di- 
mensionless curvature is k ~ 1/6.4. 

The evolution function of stiffness modulus and compres- 
sive yield strength can be get by the results of Suiker [? ] 
(a)The function of stiffness modulus Ê, (Pa) varies with cur- 
ing time t(s). Where the R?—value of this approximation is 
0.995. 


Ê, (t) = 78100 + 20t, (39) 


(b)The function of compressive yield strength G,,,.(Pa) varies 
with curing time t(s). Where the R?—value of this approxi- 
mation is 0.991. 


ps (t) = 5984 + 2.45. (40) 


Under the influence of the parameters of printing processes 
and geometric, the buckling length of the rectangular wall 
is 18. = 0.210m when simply supported on the side, and 
lf, = 0.431m when full-clamped on the side, the plastic 
collapse length is l, = 0.421m. Therefore, the actual fail- 
ure size range of the rectangular wall should be 0.210m < 
ler < 0.421m. Further, the rotational stiffness provided by 
the side wall to the main wall is determined by the aspect ra- 
tio b/d = 2.5, so as to obtain a more precise buckling length 
ler = 0.269m. 

The dimensionless buckling length of cylindrical wall is 
ler = 13.3, which is ler = 0.578m, under the condition 
of identical parameters. According to Eqs.(31) and (28) can 
get dimensionless length of the plastic collapse lp = 1.393, 


which is lp = 0.421, when E = 0.282. That is, when curva- 
ture k ~ 1/6.4, the cylindrical wall will cause plastic collapse 
due to reaching compressive yield strength. 

When the buckling length be the same, i.e., ler = 0.269m 
and ler = 6.184, substituting ler into Eq.(43) and can get the 
curvature is k = 1/17.74. In other words, the cylindrical wall 
can have a larger frame than the rectangular wall when it rep- 
resents the same stability requirements. Obviously, compared 
with the rectangular wall, the cylindrical wall can make more 
full use of materials due to the structural integrity (equivalent 


to the full-clamped on the side) and the introduction of cur- 
vature. That is to say, the cylindrical wall can reach a higher 
number of layers and have stronger stability. 


C. Analysis of failure mechanism 


Specialy, for the given geometric parameters, when h = 
43.5mm a = 8.297, dimensionless thickness h = 1 and 
(484) % = 22.99. Therefore Ter = 22.99 - lor. 

When the stiffness modulus and compressive yield pene 
are given by Eqs.(39) and (40), £g = 
0.000256(s~1),gé = aS = 0.00041(s"'),p = 
2020(kg/m*), g = 9.81m/s?. It can be seen from Eqs.(20) 
and (31) that the dimensionless curing rate E is determined by 
the axial growth rate l of the cylindrical wall. 

Q tin h 


F h 
(= = . — = 0.00280— 41 
UnhT, 2th r ant wa i (H) 


20 
78100 


The function of dimensionless curing rate, Ep and €,, with 
dimensionless radius 7 can be obtained 


= 1 Do \ike 2 
=(= = 0.004 
fp (~ 7u T, 
z, = S2l@vol _ 0.044097, (42) 
pgi 


Substituting Eq.(42); into Eq.(25) for numerical solution, 
the dimensionless radius 7 range from 5 to 10. Where the 
R?—value of this approximation is 0.999. 


Ter = 5.63943 + 44.07928 exp(—0.285097). (43) 


The Suiker number can be calculated by Eq.(33) under such 
parameters, A = 6.94196. The function of plastic collapse 
length lp with dimensionless radius 7 can be obtained by 


combining Eqs.(28) and (42)2. Then plot Î, (7)-A and Eq.(43) 


in Fig.14. 
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Figure 14: The plastic collapse length lp and elastic buckling length 
ler of 3D printed cylindrical wall vary with the dimensionless radius 
T under the given parameters. 
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When the printing parameters 
Eqs.(39)(40) and Table I, 


fixed, see Fig.14. The intersection point A of 1, (7) - A and 


are determined by 


the Suiker number A is also 


ler(T) represents the critical point of the competitive rela- 
tionship between plastic collapse and elastic buckling, i.e., 
ler = lp. At the intersection point, the dimensionless radius 
is T = 7.73496, the failure length is ly = 10.53515, which is 
equivalent to the curvature is k = 1/7.73496, and the failure 
length is lp = 0.4582m. When the dimensionless radius 
T < 7.938, the plastic collapse is the failure mechanism, 
otherwise the elastic buckling. 
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Figure 15: The area division diagram of effective length and failure 
length of 3D printed cylindrical wall under the given parameters. 


The zones of the effective length, failure length, elastic 
buckling and plastic collapse can be distinguished by Suiker 
number, see Fig.15. Blue region B is the effective length, 
yellow region C is the plastic collapse, red region D is the 
elastic buckling. In the printing processes of cylindrical wall, 
the expected layer can be achieved by selecting the curvature 
or adopting the corresponding auxiliary means to enhance the 
structure in different failure areas. 


VI. CONCLUSION 


In order to improve the overall stability of the 3D printing 
straight wall, a simple and feasible method is proposed in this 
paper, that is, the curvature along the direction of the contour 
is appropriately increased to print wavy wall. The shape of 
single wave can be cylindrical or not, and the curvature is vari- 
able. This paper based on the work of Suiker et al., extends 
the straight wall to wavy wall by increasing the curvature, and 
makes a detailed study of stability of cylindrical wall to ver- 
ify this method. It is found that the critical buckling length of 
cylindrical wall is more than twice that of the rectangular wall, 
indicating that increasing the curvature of straight wall struc- 
ture along the contour direction is a reasonable solution to im- 
prove the overall stability. From mechanics point of view, this 


method essentially improves the overall stability by introduc- 
ing curvature to enhance the overall topological stiffness of 
the straight wall, and wavy wall is a simple and feasible con- 
figuration. Specifically, the following five conclusions can be 
summarized: 


e The parameter model essentially describes the varia- 
tion characteristics of buckling length of heterogeneous 
cylindrical wall under subjected to self-weight in the 
printing processes. 


e In the specific study, each wave was treated as a cell. 
When the cell was taken as a cylindrical wall, it was 
found that the buckling length was highly sensitive to 
curvature k. With the decrease of curvature, the same 
increment of dimensionless curing rate dE will lead to 
greater increment of buckling length dler. 


e Fig.9 and Fig.10 describes the change trend of length 
scale ratio curve with curing rate, which can be used 
to study the competitive relationship between two fail- 
ure mechanisms. The Suiker number A determined by 
Eq.(33)3 describes the boundary between plastic failure 
and elastic buckling. When 1,,./ Ls > la > A, is plastic 
failure, and versa is elastic buckling. The curing ratio 8 
has a great influence on the change trend of this curve. 
Where A is determined by the initial parameters, and 
the stiffness curing ratio 8 is determined by the curing 
parameters. These two parameters determine the domi- 
nant position of elastic buckling and plastic failure. 
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